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Abstract 

We consider an atomic beam reservoir as a source of quantum noise. 
The atoms are modelled as two-state systems and interact one-at-a-time 
with the system. The Floquet operators are described in terms of the 
Fermionic creation, annihilation and number operators associated with 
the two-state atom. In the limit where the time between interactions goes 
to zero and the interaction is suitably scaled, we show that we may ob- 
tain a causal (that is, adapted) quantum stochastic differential equation 
of Hudson-Parthasarathy type, driven by creation, annihilation and con- 
servation processes. The effect of the Floquet operators in the continuous 
limit is exactly captured by the Holevo ordered form for the stochastic 
evolution. 

1 Introduction 

Periodically kicked quantum systems are described dynamically by applying a 
unitary V, called the Floquet operator, every t seconds. In an open systems 
model, V will be a unitary on a Hilbert space [}s (8) ijR where l)s is the system's 
state space and f)i? is the state space for the environment at that particular time. 
Averaging over the environment leads to a dissipative reduced dynamics on ijs- 
Here we address the question of an system being periodically kicked by inde- 
pendent environments and consider the continuous time limit r ^ 0. To obtain 
an limit open dynamics we must re-scale the Floquet operators appropriately. 

The approach which best captures the limit dynamics is that of the "time- 
ordered exponentials" introduced by Holevo pp. This theory essentially deals 
with quantum stochastic Floquet operators and is equivalent to the usual Hudson- 
Parthasarathy approach 2 . We prefer the terminology Holevo ordered from for 
the former and Wick ordered form for the latter. This is because time ordered 
exponentials appear as Dyson series expansions: here the chronologically or- 
dered terms lead to a Weyl-Stratonovich theory as opposed to the Wick-Ito of 
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Hudson and Parthasarathy [3] . Previously, we have established a quantum cen- 
tral limit for time-ordered exponentials involving reservoir fields (A) , (A) 
satisfying commutation relations [a~ (A) , (A)] = ^G'(^^), where A was 
a small parameter, and emission, absorption and scattering where present 0]. 
This is interpreted as Markovian limit where the auto-correlation time r cc A^ 
vanishes. The fields af (A) are actually superpositions of creation/annihilation 
fields for fixed momenta states modulated by a t-dependent phase component: 
the Markovian approximation is then an infinite bandwidth limit. 

In the present situation, the fields are in discrete time and are Kronecker- 
delta correlated. In this sense they are already discrete white noises. However 
we use the same strategy of anticipating the limit in terms of suitably scaled 
collective operators leading to a quantum central limit: see and chapter II 

of [7] . We have the advantage here that the dynamical updates involve operators 
independent of past state of the reservoir and so we avoid the finite-memory 
features appearing in Markovian approximations. 

The limit we consider has been used to describe the open dynamics of a 
laser mode interacting with an atomic beam reservoir |B] . This problem has also 
been studied recently by Attal and Pautrat and we obtain similar results to 
theirs. They, however, investigate the vacuum limit and use the Guichardet's 
representation of Fock space processes and the toy-Fock approximation to Fock 
space [7j. For our purposes, we find the connection to Holevo's formahsm the 
most transparent. We conclude with a construction of non- vacuum limits. 

2 Model For an Atomic Beam Reservoir 
2.1 Open Floquet Dynamics 

Let be a unitary operator on the Hilbert space [)s ^ f)_R where f)5 is states 
space for a system of interest and the state space for its current environ- 
ment. If we fix a reference density operator for the environment, then a 
completely positive map, S. on the algebra B {i)s) of bounded system operators 
is determined by 

tv^, {qs S (X)} = tvt,,^^^ {gs<^QR [X (g, 1) V} (1) 

If the environment was ignored, then V would be referred to as a Floquet 
operator, particularly, when applied repeatedly. In such cases S (X) = XV 
is a closed, and so non-dissipative, evolution. 

Our aim is to study open Floquet systems and for simplicity we shall assume 
that the repeated applications of the Floquet operator involve copies of the same 
unitary however with different, independent environments. 

We consider a repeated interaction strategy as a discrete-time open dynam- 
ics. At times t — t, 2t, 3t, . . . we have an application of a copy of the Floquet 
operator V. Let \)ii^k be the state space at time t = kr - this will be a copy of 
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- then we are interested in the Hilbert spaces 

[t/r] oo 

^l=<S)^R^k, -5^= (g) i)R,k, ^r = ^1®^'r (2) 

fc=l fc=[t/Tj+T 

where [x\ means the integer part of x. (We fix a vector eo G f)fi and use this to 
stabilize the infinite direct product.) We refer to and Sj^^ as the past and 
future; reservoir spaces respectively. 

The Floquet operator to be applied at time t = kr will be denoted Vk and 
acts on the joint space f)s but has non-trivial action only on the factors f)s 
and f)ij,fe. The unitary operator U^^^ describing the evolution from initial time 
to time t is therefore 

U^'^ = V^tM---V,Vi (3) 

It acts on 1)5 ® i^Jj but, of course, has trivial action on the future reservoir 
space. The same is true of the discrete time dynamical evolution of observables 
X €B (f)s) given by 

4^\X) = UI^^HX^I)UI'\ (4) 
We therefore have the difference equation 

7 (f^&lj+r - ^wri) = (nv.j+r - 1) u[:ii ■ (5) 

Our objective is to obtain a (quantum) stochastic differential equation for lim- 
iting situation r ^ 0. To this end, we require a reference state for the reservoir 
and we choose the pur state determined by the vector ^'^ on given by 

= eo (?) eo ig) eo (g) eo • • • 

and, since eo will typically be identified as the ground state on we shall call 
the vacuum vector for the reservoir. 

The situation can be describe alternatively as follows. The Hamiltonian 
describing the combined system and reservoir is the formal operator oni)s 
given by 

Ht^=f26{t-kT)Hl:^ (6) 
fc=i 

where h')^^ acts non-trivially only on the factors t)s and \)R^k- The Floquet 
operators are then 

14=exp{-iTWi")}. (7) 

Note that we include a dependence on the time-scale parameter r in as we 
shall require some control over the interaction in the limit r — > 0. 
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2.2 Atomic Beam Reservoirs 



One situation that wc can model in this way is when the reservoir i? is a beam 
consisting of discrete atoms, each having state space f)/{. The system might be 
a photon mode inside a cavity. The atoms pass through the cavity in a regular 
sequence and its is assumed that the atom-mode interaction takes place over a 
time period Tint shorter than the time r taken for a single atom to pass through 
the cavity: therefore the atoms are independent and at any time at most one 
atom interacts with the mode. We shall therefore assume that the interaction 
is instantaneous - that is, the system receives a "kick" from each atom. 

For simplicity, the atoms are taken to be just two-level atoms with ground 
state cq and excited state ei. The transition operator from the ground state to 
the excited state of the k—th atom is and is the ampliation of a"*" = |ei) (eo|: 
its adjoint is denoted as a~ = |eo) (ei| and the ampliation by a^. The operators 
(7^ are Fermionic variables and satisfy the anti-commutation relation 

Wt,4} = o^ Wk^4} = '^- (8) 

The operators commute for different atoms. The Hamiltonian for the beam is 
(formally) Hr = Efeli ^CTfeC^fe ■ 

The preparation procedure is the same for each atom and corresponds to an 
ensemble state g on We take the general form 

g — Picr^ + Pocr^ (9) 

where pi,pa are the probabilities to be in the excited state and ground state 
respectively. The vacuum state is, of course, specified by po = 1. 
The Hamiltonian is specified by setting 

Hk := -Hii O a-^a^ + -^Hw (8) + -^Hm ®cFk+ ^oo- (10) 

where we take Hu and i?oo to be self-adjoint and require that (i?oi)^ = -f^io- 

Wc may identify Hqq with the free system Hamiltonian Hs, while Hn may 
be considered to contain the Hr as a component. We shall assume that the 
operators are bounded with Hu also bounded away from zero. 

We shall also employ the following summation convention: whenever a re- 
peated raised and lowered Greek index appears we sum the index over the values 
zero and one. With this convention, 



were we interpret the raised index as a power: that is, [a;]° = 1 





(11) 
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[A- {t;T),A+ (s;t)] 



2tA {t As,t), 



2.3 The Collective Operators 

We define the collective operators A^ {t; r) , A {t; r) to be 

[t/rl lt/r\ 

A+ (t; t):=v^y1 '^fe ; ^" ^) •= E ; 
fe=i fe=i 

[t/r] 

A(t;r):= ^ a+a^. (12) 

fc=i 

For times t, s > 0, we have the commutation relations 

'tAs 

T 

[A{t;T),A+{s;T)] = A+{tAs;r), 
[A-{t;T),A{s;T)] = A-{tAs;T), 

where s At denotes the minimum of s and t. In the limit where r goes to zero 
while s and t are held fixed, we have the approximation 

[A- {t; t) , A+ (s; r)] w i A s. (13) 

This suggest that the collective fields A^ (t; r) converge to Bosonic quantum 
Brownian motions as r — > and that A {t; r) will converge to the Bosonic 
conservation process. 

2.4 Bosonic Noise 

Let f) be a fixed Hilbert space. The n-particle Bose states take the basic form 
(j)^® ■ ■ ■ 'S>4>n = X^crge 4'i® ' ' ' ® 4'n where wc sum over the permutation group 
&n- The n-particle state space is denoted f)®" and the Bose Fock space, with one 
particle space f), is then r+ (1^) := ^i*'^ vacuum space f)®'^ spanned 

by a single vector ^. 

The Bosonic creator, annihilator and differential second quantization fields 
are, respectively, the following operators on Fock space 

(-0) '?!>i«) • • • = Vn + lip(g)(l)j^(E) - ■ ■(S)<j>„ 
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where ^ e\) and T eB ([)). 

Now choose \) = (M+ , dt) and set 

^± := A± (l[o,t]) ; At := (i[o,t]) (14) 
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where lio,t] is the characteristic function for the interval [0,t] and l[o.t] is the 
operator on (K"'', dt) corresponding to multipUcation by l[o,t]- 

An integral calculus can be built up around the processes Af , At and t and 
is known as (Bosonic) quantum stochastic calculus. This allows us to consider 
quantum stochastic integrals of the type /^^ Fio (t) (g) dAf^ + Fqi (t) (g) dA^ + 
Fii (t) (g> dkt + Foa (t) (gi dt on (jo (E> T+ (i^ (M+, dt)) where (jo is some fixed 
Hilbert space (termed the initial space). 

We note the natural isomorphism f)o <Xi r_|_ (i^ (R+,(i/;)) = l)^] (g) ()(t where 
l)t] = f)o r+ ([0,i] ,dt)) and f)(t = r+ {{t,^),dt)). A family (F*), of 
operators on f)o (g) r+ (L^ (R^, rfi)) is said to be adapted if Ft acts trivially on 
the future space ()(( for each t. 

The Leibniz rule however breaks down for this theory since products of 
stochastic integrals must be put to Wick order before they can be re-expressed 
again as stochastic integrals. The new situation is summarized by the quantum 
Ito rule d (FG) = (dF) G + F (dG) + (dF) (dG) and the quantum Ito table 
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dA 


dA- 


dt 


dA+ 














dA 


dA+ 


dA 








dA- 


dt 


dA- 








dt 















It is convenient to denote the four basic processes as follows: 



At, 


(1,1); 


At, 


(1,0); 




(0,1); 




(0,0). 



The Ito table then simplifies to dA^^dA^'^ = except for the cases 

dAfdAl^ = dAf. (15) 

The next theorem is from 

Theorem 2.1: There exists an unique solution Ut to the quantum stochastic 
differential equation (qsde) 

dUt = L^p ® dAf, Uo = 1 

whenever the coefficients Lap are in B{i)o). The solution is automatically 
adapted and, moreover, will be unitary provided that the coefficients take the 
form 

Ln=W-l; Lw = L; = -L^W; Loo = -tH-^L^L (16) 

with W unitary and H self-adjoint. 
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2.5 Convergence of the Collective Processes 

Let (/), ijj G (R^, dt) be Riemann integrable and T € L°° (M+, dt) continuous. 
We define the pre-limit fields 

A+ i<l>,T) := ^^^^<^(rfc) a+, A" (V,r) := V^^r {rk) 

k k 

Vt/T\ 

A{T;r) := ^T(rfc) a+a^. 



fe=i 



then 



[A- (V, r) , A+ {<i>, r)] = r ^ V* {rk) 4> (rk) - 2t ^ V* (rfc) (rfc) a+a^ . 

fe k 

which converges to (V'|'/>) = ^* (i) {t) dt as r ^ in the vacuum state 
More generally we have convergence of the type 

hm f[ exp{iA(T,-,T) + a+ (.^,.,r) (V,-,t)} $M 



[| exp {iA {Tj) + iA+ + iA- (V,) } * 



(17) 



where 11^=1, •••."^j — ^1^2 ■ ■ ■ Xn denotes an ordered product of operators. 

3 Decomposition of the Floquet operators 

Let 0-^ be the two-level transition operators. We set 



n = H, 



a(3 I 



and assume that < C. Using the anti-commutation relations, we have 

that 

(r T-Cf = [HiiHii + T HiiHqo + r HiqHqi — r ifoi-ffio + -ffoo-ffii] cr+cr" 
+7? [HiiHw + tHioHoo + tHooHio] ct+ 
[ifoi^fii + tHoiHoo + tHoqHoi] u' 
+r [HqiHio + T i?oo-f^oo] • 

More generally, for m > 2, we have that (recall our summation convention!) 



(rH)" = r"i/„„^„---7J„i/3, 







a 


/3„ 


■(7+- 


ai 


a 

















(18) 
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and, again by repeated use of the anti-commutation relations, we are lead to 
the form 





a 


a 









/3 



where 



Xi;) (r) (i/n)"-" H,p + O (r) 



(19) 
(20) 



The remainder is a polynomial of degree m in r, whose coefficients are sums 
of n-fold products of the Ha/3's, having no constant term. Here O (t) means a 
term going to zero in operator-norm faster than r as t ^ 0. 
We next of all compute the Floquet operator: 



V = exp {-itH} = 1 + t {L^p + R^p (r)} < 





OL 


a 









(21) 



Here the L^fs and Raf} (t) are bounded operators on the system space with 
Raf3 (t) = O (t). Explicitly, the coefficients Lap are given by 



ill 

Lqi 
Lqq 



-iHii 



Hn 



-Hio] 



-iHi 



H, 



01" 



-iHn 



'iHi 



Hoi 



1 + iHii 



(HnY 



H. 



(22) 



We remark that these coefficients take the form (|16|l where W := exp{—iHii} 



is unitary, H Hqq — Hqi 



Hii-sin(gii) 



(Hll) 



otherwise arbitrary. (Note that — 



sm X 




-fflo is self-adjoint and L is bounded but 
> for a; > 0.) 



4 Limit For the Vacuum State 



We begin by assuming that pi = and that therefore the reservoir is in the 
vacuum state '^^ . 

(t) 

Let t > 0, then we are interested in the unitary — V[_t/T\ • • • V2V1 in the 
limit r — > 0. By virtue of an uniform estimate established in the next section, 
we have that the components Rap (t) make negligible contribution in the limit 
T ^ 0. It is easy to see that, if fc = [t/rj and if we ignore the negligible 
component , then 



1 



[Vk+i - 1) u, 



La/3 





OL 


'^k+1 
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We shall replace this by a quantum stochastic differential equation shortly. The 
operator [/j^' is then represented as 



M= n <! l+T[£a/3+i?a/3(T) 









m \ 
























\ 



(23) 



Theorem 4.1: In the above notations, the discrete time family l^t^"*! con- 
verges to quantum stochastic process \Jt on f)s r_|_ [l? (K+, rft)) in the sense 
that, for all u,v G t)s, integers n,m and for all (j^jjipj G L,^ (M+,(it) Riemann 
integrable, we have the uniform convergence 

(a+ r) • • • A+ (01, r) u ® (7f ) A+ (^„, r) • • • (Vi, r) z; ® 

{A+ (0„J • • • A+ (0i) M ® *| C/t A+ (^„) • • • (V-i) « ® *) (24) 

The process C/j is moreover unitary, adapted and satisfies the (quantum) stochas- 
tic differential equation 



dUt ^ L^f3 (E> dAf Ut, Uo = l 



(25) 



Remark 1: The solution to H25|) can be written as 



[/t = N cxp I y ds Lafi [a+] " [a^ ] 



where are quantum white noises. Here the symbol N stands for normal 
ordering and we understand the formal development 



„=l"'*>ti>-*'i>0 



dtn ■ ■ ■ dti La„f3^ ■ ■ ■ LaiP-^ 



X [at 



+ ' 



The connection with the Hudson-Parthasarathy notation is made by the 
replacements [a^]" Kt [a^]^ dt KtdA'^^ for adapted Kt- 

Remark 2: There is an alternative presentation 1 which we refer to as the 
Holevo ordered form. We write 



and understand this to be the Ito qsde 



dUt = 



-l]Uu t/o = l. 
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From the quantum Ito table we have that 



H„„dA° 



- 1 = 



E 



H. 



. Ha,f3^ dA 



.dA. 



ai/3l 



Ha/s dAf 



Y.-^H^i{H,,r-^H,p dAf. 



We see that the Ito coefficients, Lap, and Holevo coefficients, Hap, are connected 
according to the same relations as (I22|) . 

Remark 3: The theorem may be restated in a more elegant fashion. The dis- 
crete unitary process Ut (r) = Texp |— i Hg (r) ds|, where Ht (r) = Yl,k ^ {t — kr) Tik 
with Tik given by converges to the continuous-time unitary process Ut = 
H exp I ds Hap [atr[a-fy 

Remark 4: The basic estimates in section 4 of serve to show the convergence 
of H23|l to the Holevo ordered form. 

Next of all we turn our attention to the Heisenberg evolution. We begin by 

/9 



noting that if we set AA"^ = r 
AA^'^AA^^ = T 





a. 


(7 









then 



AA"''^ 



and that otherwise AA°''^ AA^'" = O (r) AA''?. This is the discrete form of the 
quantum Ito table (|15|) . 

Let J^^^ {X) = [/(^'■'^ {X ® 1) Ul''\ If /c Yt/T\ , then 

{4''> + (XL„;3) -I- ) (iLXii^) + (r)} ( 



'k+l 



'k+1 



where Sap (t) — O (r) . Again the terms Sap (t) will have negligible contribution 
in the r — > limit. We establish the appropriate uniform estimate in the next 
section. 

Theorem 4.2: In the above notations, the discrete time family | Jf (-^)| con- 
verges to quantum stochastic process Jt {X) — u} {X ® 1) Ut on f]5(X)r+ (L^ dtf) 
in the sense that, for all u, -y S fig, integers n, m and for all (pj^i^j G -^^ (^^, dt) 
Riemann integrable, we have the uniform convergence 

(a+ (^,„,t)...A+ (01, t) W ^+ (^„,t)...A+ (V'i,r) 

^ (A+ (0„) • • • A+ u ® vl/| J, (X) A+ ■■■A+ (Vi) « ® *) . (26) 
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The process Ut is moreover unitary, adapted and satisfies the (quantum) stochas- 
tic differential equation 



where 



dJt {X) ^ Jt iUpX) ® dA'^'^ , Jo{X)^X(g,l 
i^apX L^p^X + XLafj + L\^XLip. 



(27) 
(28) 



Remark 5 A completely positive semigroup {Sf : t > 0} is then defined on 
B (f)s) by {u (g) *| Jt (X) v(g)^) := {u\ Et (X) v) and we have Et = exp{t£oo} 
where the Lindblad generator is 

£00 (X) ^^[L\X]L + [X, L] - I [X, H] 
where L = ^'""'^ Hio and H ^ Hoo - ^01 """r^'lf ^10- 



5 Uniform Estimates 

We now want to obtain a norm-estimate for the series (|21|l based on the expan- 
sion appearing in (|18|l . 



Lemma 5.1: In the notations of the previous section 



n>0 ■ q,/3g{0,1}" 



H. 



■H, 



Pi 



< expjr [e^ - l)} ■ 

(Recall that C — maxQ,^ H^^a/sH , so e*-^ > 1, and note that we have broken from 
our summation convention to show explicitly that we have a sum of norms.) 
Proof. Evidently we have the bound 



a,/3e{0,l}" 

Q,/3e{oa}" 



■H. 



/3„ 



Pi 



(29) 



Now • 



'a+' 




a 


/5„ 




Ql 


(7 

















Pi 



can be described as follows: we have 



n ordered vertices labelled j = 1 , • • • ,n and at the fc-th vertex will be either 
a~ , ^/Ta~^ , yjTa~ or r depending on whether (aj, (i^ = (1, 1) , (1, 0) , (0, 1) 
or (0, 0) respectively. Our first objective is to put this expression to Wick order. 
In placing the ct+'s to the left of the cr^'s, we must repeatedly use the anti- 
commutation relations. This introduces, in the usual way, the notion of pair 
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contractions between a and a a~ at different vertices- that is, we replace 
Aa'^ Ba~ C with ABC . (We ignore the possible minus sign occurring as we want 
a norm estimate.) 

We then have to contend with a sum over all possible pair contractions. We 
have at most one creator and one annihilator at each vertex. Therefore, in a 
typical term we shall have several vertices connected through pair contractions 
and these vertices form disjoint subsets of all the n vertices. We also take each 
(0, 0) vertex to be a singleton set. In this way each term corresponds to a 
partition on the n vertices into subsets. Now recall that the number of ways 
to partition n objects into m subsets is given by Stirling's number, S {n, m), of 
the second kind jn|. Each subset of the partition contributes a factor r: this is 
obvious for the singletons and more generally we have a subset with ^/t for both 
of the terminal vertices and unity for the internal scattering vertices. We have 
in addition a product of the uncontracted cr* but this will have norm bounded 
by unity. 

This leads to the following bound for H29|) : J2m=i C"t™S (n, m) .The norm- 
bound for the expansion for the full series is then 

oo -. n 

E;^ E^"^"^('^'"^) = exp{r(e^-1)} (30) 

n—O m—1 

where we use the well known generating series ^01 for the Stirling numbers. ■ 
(We remark that the integer S (n, m) gives the coefficient of A"* in the n-th 
moment of a Poisson distributed random variable with intensity A.) 

Corollary 5.2: 



xn E 

j a,/3e{oa}"3 



E E 

k—l ni ,- ■ ■ ,n/^ >0 



^niH hrifc 

7ii! • • • nfc! 













Ql 















< exp{t(e'^- 1)} 



This establishes a uniform estimate for the series expansion of Ut = y\t/T\ ■ ■ ■V2V1 
based on the development (|21|l . We now do the same for the Heisenberg evolu- 
tion. For X E B {t)s) we have 



a 



01 
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and by the previous arguments we see that the sum of the norms of these 
summands can be bounded by 

r^n+n' 'm-\-m' 

ll^ll E E ^i^^M S (n, m) 5 (n', m') = exp {2r (e^ - l)} . 

n,n' m^m' 

Likewise, the series expansion of J^*'^'' (X) — U^'^^'^ {X ® 1) J/j^-* in terms of the 
fundamental Hamiltonian components wiU be bounded by ||X|j exp {2t (e*^ ~ l) } 
uniformly. 

6 Non- Vacuum State 
6.1 Gaussian Case 

We consider the situation where Hn = 0. Let (.) be the state determined by 
the density matrix g in 0. With the convention that = 2(T+c7~ — 1, we set 

= VPo^t ® 1 + VPl'^k ® cr^ 

These operators commute for different lables k and we have {o''^,^k} — 
(aj)^ — 0. Moreover, algebra generated by the with the pure state eo (8) eo 
is isomorphic to the one generated by the ctj, with mixed state (.). 

If we adopt the new representation then we can consider collective fields 

(0, r) = v^B± (0, r) + VpTC^ t) 

where 

B+(0,r) : - V^^0(fcr) 

k 

Here j : i— !■</)* is the complex conjugation. We note that 

[B- (0,r),C+ (V,r)] =2r^</>(fcr)Xfcr) a+ ® 

fe 

and that this and similar terms are negligible in the limit r — )■ for the state 

(g) 

We therefore find that the fields (.,t) and (.,t) in the state $^ (g) <I>^ 
converge in distribution to independent (that is, commuting) Bose fields i?" (.) 
and C* (.) with the double Fock vacuum state ^iSi^. We set A (.) = ^/pa B^ (.) + 
^yPlC^ (j.) and find the modified Ito table 

dA^dA^ = Po dB^rdB^ = po dt, dA^dA^^ = pi dC^ dC^ = pi dt. 
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It is relatively easy to sec that there again exists a limit process Ut, this time 
driven by the fields Af. The limit qsde will be 



Hio dAf + Hi 



01 



which is again unitary and adapted to the noise fields A: 



t ■ 



■)} 



(31) 



6.2 Asymptotic Case 

We consider a new orthonormal basis for i)R 

|eo) = 0VPo\eo) +0*VPi\ei); 
|ei) = -rV^|eo) + 0v^|ei). 

where is a complex number of unit modulus and Po + Pi = 0. Introduce new 
transition operators 

cf-+ = |ei) (eo| ,0-" = |eo) (ei| 

The new variables again satisfy the proper anti-commutation relations how- 
ever a~ annihilates the state eo- 

We note that (eo| cr+a" |eo) = pi, (eo|<T~cr+ |eo) =po, (eo| o-+ |eo) = ^/pipo^^ 
and (eo| cr~ |eo) = s/PiPoO*"^- We find that 

(7+ = - pid~ - y/pipo {0^ + 9*^) a+a~ + ^/plPo&^^, 

a+a- = ipo-pi)a+a- -./pm{0'^ + 0*'^){a+ + a-)+pi. 

After some rearrangement, we obtain 



H, 



a/3 



1/5 



= Ha0 O 



a r ~ - 

a 



where 
Hoi 



= H, 



00 



PaPl a2 



f iioi H -"11, 

T T 



= poHiQ - piHqi + 



= poHoi - piH 
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POPl „*2 TT 

f iiii; 



{po - Pi) - vwT (r + 9*^) Hoi - VwT (r + 9*^) H 



Let us now adopt cq as the new ground state and stabilizing vector. We set 
= eo (8) eo (8) eo (8) • • • . Clearly we would like to use the tilded variables in place 
of the un-tilded ones however there is an explosion problem associated with the 
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coefficients Hafj. One way to resolve tfiis is to allow the state to depend on r 
by taking 

pi =7^ + (t2) . (32) 
In this case we have the asymptotic behaviour 

ffoo - i/oo + T^^'^io + 7^*'^oi + + O (r) ; 

Hio = i/io+7^'^ii+0(r); 

Hoi = Hoi+je*^Hn+0{T); 

Hn = Hii + 0{t). 

Here the continuous limit follows by using the above forms for the Haf3 and 
ignoring the O (r) terms. The scaling used in (|32|l is necessary if we wish to 
obtain a Gaussian limit for the collective operators: this is related to the notion 
of macroscopic states in statistical mechanics 
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